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The quantum scattering of magnon bound states in the anisotropic Heisenberg spin chain is shown 
to display features similar to the scattering of solitons in classical exactly solvable models. Local¬ 
ized colliding Gaussian wave packets of bound magnons are constructed from string solutions of 
the Bethe equations and subsequently evolved in time, relying on an algebraic Bethe ansatz based 
framework for the computation of local expectation values in real space-time. The local magneti¬ 
zation profile shows the trajectories of colliding wave packets of bound magnons, which obtain a 
spatial displacement upon scattering. Analytic predictions on the displacements for various values 
of anisotropy and string lengths are derived from scattering theory and Bethe ansatz phase shifts, 
matching time evolution fits on the displacements. The time evolved block decimation (TEBD) 
algorithm allows for the study of scattering displacements from spin-block states, showing similar 
scattering displacement features. 


I. INTRODUCTION 

The study of classical dynamics in nonlinear media has 
proven to be a source of astonishing surprises over the last 
century. Two observations, both based on numerical sim¬ 
ulations, have challenged prejudices and fundamentally 
altered traditonal ways of thinking. First, the famous ob¬ 
servation of Fermi, Pasta, Ulam, and TsingorPof a simple 
nonlinearly-coupled set of oscillators showing nontrivial 
recurrences, has shattered the long-held assumption that 
all nonlinear dynamical systems ergodically explore their 
full phase space. Second, thepioneering numerical anal¬ 
ysis of Zabusky and KruskaP on the Korteweg-de Vries 
equatiorP demonstrated that this equation supports ex¬ 
citations, which they coined ‘solitons’, displaying a num¬ 
ber of surprising fundamental features. The solitons are 
localized in space, with a form remaining stable under 
time evolution which sees them moving uniformly at a 
speed linearly proportional to their amplitude. Addition¬ 
ally the astounding characteristic was observed of solitons 
emerging intact from mutual scattering processes, during 
which they simply ‘ “pass through” one other without los¬ 
ing their identity’,1^ the only effect of the collision being 
a relative spatial displacement as compared to their free 
propagation. The proper understanding of solitons in 
nonlinear media ultimately led to t he d evelopment of the 
classical inverse scattering method,®^ which is the over¬ 
arching framework for classical integrable models. 

Nonlinear classical systems find their quantum me¬ 
chanical analogue in the shape of interacting many-body 
systems. In parallel to the classical case, some quan¬ 
tum models have been shown to be special, in the sense 
of being exactly solvable using the quantum version of 
the inverse scattering methocP (one might also say in¬ 
tegrable, though the quantum notion of integrability is 
not as well-defined as its classical counterpart^) . Funda¬ 
mental representatives of this family are the Heisenberg 


spin chain, solved by Bethe using what is now known as 
the Bethe ansatz,along with the Lieb-Liniger model of 
(5-interacting bosons on a line.l^ For the latter, particle¬ 
like excitations called Lieb Type I modeJi^ exist due to 
the interparticle repulsion, along with Type II hole-like 
excitations visualized as holes in an effective Fermi sea. 
It is possible to distinguish the presence of Type I and 
II modes in correlated b osonic gases in optical lattices 
using Bragg spectroscopy.l^^J^ Quantum magnets such 
as the Heisenberg spin chain similarly carry particle-like 
magnon modes when the magnetization is close to satura¬ 
tion. In the limit of small magnetization, hole-like modes 
again appear, which in zero held are known as spinons.^ 
Their dynamics can be e xperim entally observed using in¬ 
elastic neutron scatteringThe Heisenberg chain sup¬ 
ports distinct bound states of magnons, whose dynam¬ 
ics has been investigated theoreticallj^^ and has recently 
been observed experimentally.l^^ 

One could view such excitations as the quantum equiv¬ 
alents of classical solitons. This equivalence is however 
only partial: on the one hand, these quantum mechanical 
modes represent exact eigenstates and are stable under 
time evolution; on the other hand, being exact eigen¬ 
states of translationally-invariant systems, they are not 
spatially localized. That said, as is usually the case in 
quantum mechanics, it is possible to adopt a ‘complemen¬ 
tary’ picture and create spatially localized excitations by 
forming wave packets of fundamental excitations by lin¬ 
early combining states over a range of differing momenta. 
Locality however comes at a price: the wave packet, mix¬ 
ing together states at different energies, will disperse and 
is thus not stable over long timescales, unlike its classi¬ 
cal counterpart. We use the term quasisoliton for such 
a wave packet construction, an example of which was re¬ 
cently studied in the context of the Lieb-Liniger model,^ 
while their mut ual sc attering has been studied in quan¬ 
tum spin chains 
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The spectroscopic methods traditionally employed to 
experimentally study condensed-matter systems typi¬ 
cally provide momentum- and energy-resolved measure¬ 
ments. However, current experimental developments pro¬ 
vide motivation to obtain a better theoretical under¬ 
standing of spatially localized dynamics. Time-resolved 
experiments are now able to track quantum many-body 
systems at timescales smaller than relaxation tim escales, 
particularly in experiments with ultracold atoms,^l ^^ l ^^ l 
but also in NMR setups,!^ and potentially in pump-probe 
spectroscopy experiments.!^ As a result, non-dissipative 
dynamics of many-body systems out of equilibrium is 
now a rapidly growing field of experimental as well as 
theoretical investigation.!^ Novel cold-atom experimen¬ 
tal techniques for spatially resol ved man ipulation and ob¬ 
servation at the single-site leve P^ ! ^^ -*^ have opened the 
door to explicit high-resolution tracking of spatial prop¬ 
agation phenomena. Moreover, an experiment with in¬ 
teracting bosonic atoms has highlighted the interaction- 
induced longevity of repulsive pairs.!^ This has moti¬ 
vated increased theoretical attention to the (anti-)bind- 
ing of localized excitations and inter action s of these 
bound clusters, b oth in itinerant systemJiSlll! and in spin 
chains .!i^!i2lM29l Very recently, the spatial dynamics of 
itinerant clusters has been studied in experiment.!^ Prop¬ 
agation of quantum solitons in Bose-Hubbard chains has 
also been studied numerically.!^ 

In Ref. [ini the scattering of a magnon wave packet on 
approximate bound eigenstates of n particles was studied 
numerically in the Heisenberg chain. In the present work, 
we present an algebraic framework and exact calculations 
based on Bethe ansatz. We therefore consider quantum 
scattering of localized excitations over a ferromagnetic 
background in the Bethe ansatz sol vable an isotropic spin- 
1/2 Heisenberg chain (XXZ model)^ ^^ ! ^^ ! 


H = J 


N r 

E 




QX QX _\_ CV qV _l a I 



( 1 ) 

The XXZ model is experimentally realizable, for example 
in the setups of Refs. [T7| and which use two hyper- 
fine states of bosonic atoms in the Mott phase to ex¬ 
perimentally realize the spin-up and spin-down states. 
The effective model is a nearly isotropic (Awl) Heisen¬ 
berg chain, while an experimental setup with variable 
anisotropy A is under development.!^ In addition, the 
XXZ model has been shown to describe Josephson junc¬ 
tion arrays of the flux qubit type,!^ and may also be 
realizable in optical lattice^^ or with polaritons in cou¬ 
pled arrays of cavities.!^ It is conceivable that these or 
similar experimental setups may provide time-resolved 
observations of propagating and interacting localized ex¬ 
citations. 

The parameter J in Hamiltonian Q is given by the 
exchange interaction of two neighboring electrons or, in 
the aforementioned experimental setup, by the exchange 
interaction between two neighboring bosons in an opti¬ 
cal lattice at unit filling. Its sign does not matter for 


our purposes hereI1212S1291 for definiteness, we set it to 
J > 0. We distinguish two regions for the anisotropy 
parameter A, namely the planar xy (|A| < 1) or axial 0 
(|A| > 1) cases. Due to the interaction term (oc A) the 
XXZ chain displays a whole zoology of fundamental ex¬ 
citations: isolated down spins can form spatially bound 
states (as was understood by Bethe already in his original 
publicatiorP) , whose bond size decreases as A increases. 
These excitations are often referred to as ‘string states’, 
as the rapidities describing Bethe states containing such 
bound multi-magnons appear as approximately equally- 
spaced vertical strings in the complex plane (the precise 
set of available bound states depends on the value of A; 
at the isotropic point A = 1 and in the axial regime, all 
string lengths are allowed). 

We construct spatially localized wave packets of n 
bound magnons using linear combinations of these string 
states of length n with Gaussian-distributed momenta. 
We call them ‘n-string wave packets’. We further 
demonstrate that exact methods based on the algebraic 
Bethe ansatprovide a framework to evaluate the time- 
dependent expectation value of the local magnetization 
{Sj{t)) algebraically, which can be used to track those 
localized magnon-like wave packets. We investigate their 
stability and mutual scattering using a combination of 
scattering theory, Bethe ansatz, and numerically exact 
calculations. 

At large anisotropy A 3> 1, magnon bound states re¬ 
semble having downturned spins on neighboring sites. 
An n-string wave packet is thus closely approximated 
by a consecutive block of n downturned spins. While 
this correspondence breaks down at smaller A, this pro¬ 
vides motivation to study the evolution of states with 
downturned spins on a consecutive block of sites. In ad¬ 
dition, this is exactl y the type of initial state prepared 
in experiments! ^^ ! ^^ ! Scattering of such blocks has been 
explored in Ref. Ea where a spatial displacement of two 
sites was observed for the block at several A, and ex¬ 
plained at large A in terms of energy conservation. In the 
present work, we connect scattering phase shifts with tra¬ 
jectory displacements in order to provide a Bethe ansatz 
derivation of the observed displacements. 

The paper is organized as follows. In Sec. we will 
introduce the concepts of string solutions and the scat¬ 
tering phase shifts associated with n-strings, along with 
details on the algebraic Bethe ansatJ^ based evaluation 
of the time-dependent expectation value of the local mag¬ 
netization {Sj{t)), which can be used to track localized 
magnon-like excitations of the spin chain. In Sec. m we 
elaborate on the construction, stability and time evolu¬ 
tion of scattering of n-string wave packets. We also com¬ 
pare with the stability of consecutive-site spin blocks. In 
Sec. |IV] we consider the scattering trajectory displace¬ 
ments, derive analytical results for them and compare 
with numerical measurements. The A —>■ oo limit is 
treated analytically, comparing with numerical data for 
spin block initial states. The appendices provide details 
involving scattering theory and details on obtaining the 
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phase shift directly from the phase of the time-evolving 
wave function. 


II. TIME EVOLUTION FROM BETHE ANSATZ 


In this section we present the basic formulas of the 
Bethe ansatz for the spin-1/2 XXZ model and explain 
the concept of string solutions and the scattering phase 
shift associated with n-strings. In Subsec. II C[ expres¬ 
sions are given for the time evolved expectation value of 
the local magnetization by using results from algebraic 
Bethe ansatz. The method used here relies on the avail¬ 
ability of determinant expressions for matrix elements 
between Bethe states.SSHUl This last subsection is rela¬ 
tively technical and could be skipped on first reading. 


A. Coordinate Bethe ansatz for the XXZ model 


The eigenstates of the XXZ spin chain Q can be con¬ 
structed via the Bethe ansatz pES and have the form 


M 

iw)= E 

jl<...<jM Q a=l 

( 2 ) 

where M denotes the number of downturned spins and 
therefore fixes the magnetization. The sum over Q is 
a sum over all permutations of M objects and the am¬ 
plitudes Aq are related to the scattering phases. The 
set of M complex rapidities {A} = completely 

determines the Bethe state and is simply related to the 
physical energy and momentum. By imposing periodic 
boundary conditions, i.e. = Sf for a = x,y,z 

in the Hamiltonian 0: each set of rapidities {A} cor¬ 
responding to an M-magnon eigenstate must obey Bethe 
equations, 

The different definitions of (j)n{X), dn(X) and for various 
regions in anisotropy A are given in Tab.[l] Within Bethe 
ansatz, the momenta of single downturned spins can be 
parametrized in terms of rapidities, 


p(X) = —i In 


</i(A) 


= TT - 0i(A). 


( 4 ) 


By invoking Schrbdinger’s equation iJ|{A}) = i5|{A}) 
for a Bethe state consisting of a single downturned spin, 
the magnon dispersion relation 


E{p) = J{cos{p) - A) (5) 


is easily derived. In the case of just two single magnons, 
their scattering phase shift x can be obtained from the 


4>n{X) 


e„(A) 


|A| <1 

acos (A) sinh (A -|- 

2 atan 

f tanh(A) 

Vtan(^^) 




/2A\ 

A = 1 

A+f 

2 atan — 
\nj 

A > 1 

acosh(A) sin(A-|-^) 

2 atan | 

^ tan(A) 

ytanh(^) 


TABLE I. Definitions of functions appearing in Bethe ansatz 
for different regions in anisotropy A. 


permutation of two magnons in the Schrodinger equation, 

Aq, i + ehpi+P2) _2Ae»P2 

Aq 1 + ehpi+P2) _ 2Ae*Pi ^ ’ 

where Q is the identity and Q' interchanges the two 
indices 1 and 2. Furthermore, the magnon momenta 
Pi = p(Ai) and p 2 = p{X 2 ) as well as the scattering phase 
shift X = ^* 2 (Ai — A 2 ) are parametrized by the two rapidi¬ 
ties Ai and A 2 . 


B. Strings and magnons 


The sets of rapidities solving Bethe Eqs 0 are self¬ 
conjugate and arrange themselves in patterns of string 
solutions, 

Ea “ + 1 - 2a) -I- f-(l - z/j) -I- , (7) 

where the string center A^"^ S K and a = I,..., n is the 
internal label of a rapidity within a string of length n and 
parity vj. In the planar regime | A| < 1, periodicity of the 
trigonometric functions also allows for string centers to 
be located on the line i7r/2, resulting in negative parity 
strings (zy, = —I). This type of strings will be left out 
of consideration for the analysis of scattering magnons, 
restricting to Vj = 1. 

At finite size, solutions are not exactly given by strings, 
but rather contain string deformations (5^”^ € C, under 
the constraint that the full set of rapidities {A} remains 
self-conjugate. In the cases considered here the devi¬ 
ations are exponentially small in system size and it is 
therefore sufficient to take the limit of vanishing devia¬ 
tions. In this limit, the prodnct of the Bethe equations 
of all rapidities within a string reduces to the Bethe- 
Gaudin-Takahashi equations,!^ which are similar to the 
Bethe equations but given in terms of the n-string cen¬ 
ters A^"^. In logarithmic form they read 


1 

^ri{Xj ) A^, )- I- , 


m k=l 


/ = 1,..., M„ 


( 8 ) 
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where denotes the number of n-strings present, sat¬ 
isfying The logarithmic scattering ker¬ 

nels 9n{\) are defined in Tab. [I] and the scattering phase 
shift between two individual strings of arbitrary length 
is given by (for |A| < 1 we consider only strings with 
positive parity) 

0nm('^) ~ (1 ^nm)9\n—m\{^) T ‘^^\n—m\+2i^) 

-|- . . . -f 20„+m-2(A) + 9n+mi^) ■ (9) 

The logarithmic form of the Bethe-Gaudin-Takahashi 
equations allows for the introduction of string quantum 
numbers ^, obeying an exclusion principle for all Bethe 
states, meaning that every Bethe state is characterized by 
a unique set of string quantum numbers. By considering 
the limit of sending a string center to infinity, the max¬ 
imum allowed string quantum number can be derived. 
These limiting quantum numbers define the dimensions 
of sub-sectors of the Hilbert space containing a specific 
string content. 

Moreover, in the planar case |A| < 1, the existence of 
strings with a specific length n in the spectrum is deter¬ 
mined by the anisotropy.*^ Therefore restrictions on the 
availability of n-string wave packets as well as on their 
momentcP^ are present in the planar case |A| < 1, de¬ 
pending on the value of A. 

We get all Bethe states with the desired string content 
by solving the Bethe-Gaudin-Takahashi equations Q us¬ 
ing an iterative algorithm for all combinations of allowed 
string quantum numbers. After obtaining the rapidities, 
the energy of a Bethe state containing strings is easily 
computed as 

i?{A}=-2l'/>2(0)|^^C(A5”^), (10) 

n j=l 

where 9'^ is the derivative of The energy contribution 
of a string of length n to the energy E^x} is 

= j|^(cos(p)-e„). (11) 

02n(O) 

where e„ = cos(nC), l,cosh(n^) for |A| < 1, A = 1, and 
A > 1, respectively. The momentum p = p("i(A) of an 
n-string with string center A is given by 

p(")(A)=^-0„(A). (12) 

In the following we use the convention —tt < < tt. 

Note that for a single n-string the minimum of the energy 
dispersion is always at A*-"i = 0, i.e. at momentum = 

TT. The total momentum of a Bethe state can be extracted 
from its string quantum numbers, 

n n,j 


string quantum number can be associated to a single 
particle momentum pj”^ 


= p(")(A"^)=7r- 


^ r(n) 

N ^ 


(14) 


To summarise, the rapiditites belonging to each eigen¬ 
state are obtained by iteratively solving the Bethe- 
Gaudin-Takahashi equations Q. They can be used to 
evaluate the determinant expressions (18) for the nor¬ 
malised matrix elements of the following section. 


C. Magnetization expectation value 


Time evolution of the expectation value of the local 
magnetization {Sj{t)) is performed by making use of the 
algebraic Bethe ansatz.l^The time dependent wave func¬ 
tion is computed using unitary time evolution in a basis 
of Bethe states |{A}), 

|4/(t)) =e-^‘|4/(0)) =^e-*^tM*q,}|{A}), (15) 

{A} 


where the coefficients = ({A}|4'(0)) are determined 


by the initial state |4'(0)), which is given in Sec. Ill 


Eq. (23), for the construction of n-string wave packets. 

The expectation value of the local magnetization at 
site j is given by 


isnt))= E 

{A}.M 

(16) 

As the states |{A}), |{/r}) are Bethe states determined 
respectively by sets of rapidities {^j}pLiJ that 

obey Bethe Eqs ([^ (here we have M rapidities in both 
sets since the operator Sj does not change the magneti¬ 
zation), the matrix elements are given by the normalised 
expressions obtained from algebraic Bethe ansatz 




vmiAwm)' 


(17) 


Here we make use of the determinant expressions ob¬ 
tained in Ref. 02 




TT 4>li.9'k) 
7’i-i({A}) (/)i(Afc) 


det [g(M,{A})-2P(M,{A})] 

M 

n — Mi)<(’o(A/— Afc) 

k,l=l 

k<l 


(18) 


In the thermodynamic limit and for a finite number M of with Pj({A}) = e and 4>n{F) defined in Tab.^ The 

rapidities, each separate bound magnon represented by a entries of the matrices H{{p\, {A}) and P({p}, {A}) are 
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given by 


{A}) = - Xb) 

k^a 
N M 

0_2(/ife - Ah) , 

k^l 
k^a 

(19) 

Pabm,{X}) = —^^^^^Mh(.X,-Xt). ( 20 ) 

The normali zation A/'({A}) is computed from the Gaudin 

determinant ESESl 

A^({A}) = [</>2(0)r fi — ■ ^'~^| det$({A}), (21) 

k^l 

where the Gaudin matrix is given by the Jacobian of the 
Bethe equations, 


<^-i(Ab) 

MXb) 


M 


$ab({ A}) = Sab N0[ (Xa) - ^ 0^(A, - Afc) 


fe=l 


+ ^2(Aa ~ Ah) . 


( 22 ) 


The time-dependent expectation value of the local 
magnetization is then obtained by evaluating the double 
sum over matrix elements in Eq. (16). By construction, 


the double sum only includes eigenstates with the same 
particle content, which is not large for the few-magnon 
states we will consider. As a result, the double summa¬ 
tion is still tractable at lattice sites N ^ 0(10^). In 
the case of dealing with string solutions for the magnon 
bound states, reduced determinants for strings described 
in Ref. iTTl must be used. 


III. BOUND MAGNON WAVE PACKETS 


The strings described in the previous section do not 
correspond one to one to localized bound states of down- 
turned spins, but rather are translationally invariant con¬ 
stituents of Bethe eigenstates. In order to create states of 
n bound magnons with localized magnetization features, 
we construct Gaussian wave packets by summing over 
single n-string states (labeled by the string center A^")) 
with momenta distributed around p, 

1^(0)) = Afo^e-*P®-T(p-p)"|A(-)(p)), (23) 

p 


where Afo is a normalization constant. Unitary time evo¬ 
lution under Hamiltonian 0 implies an expression of 
the velocities of the wave packets by expanding the dis¬ 
persion relation (11) around p = p to first order (see 


Appendix [b|). 


9£;(”)(p) 


dp 




The prefactor can be expanded for large anisotropy. 


<(> 2 ( 0 ) 

02n(O) 


(2A)-”+i, 


(25) 


implying that wave packets constructed from higher 
strings have a lower velocity in real space. 

In the remainder of this section we will analyze the sta¬ 
bility of n-string wave packets, along with the stability of 
another form of a localized multi-magnon, a consecutive- 
site spin block. The two constructions are closely related 
to each other for large A. Subsequently, scattering pro¬ 
cesses of n-string wave packets are visualised by comput¬ 
ing the time evolution of the local magnetization. Fur¬ 
thermore, we describe a method for direct observation of 
the phase accumulated by the wave function of a finite 
spin chain during a scattering process. 


A. Stability of n-string wave packets 

The center of the wave packet x{t) and its width Ax{t) 
are respectively given by the expectation value and the 
variance of the position operator 

^ n \ 

■ (») 

j=i V 2 

In the continuum limit, the sum over all possible n-string 
momenta p = p^"^ can be approximated by an integral. 
For the width of the Gaussian wave packet in real space 
we obtain 

Ax{t = 0) := v/('k(0)|[f-x(0)]^|vI/(0)) « I. (27) 

String solutions, being associated with bound states 
of magnons with exponentially decaying wave functions, 
will add exponential terms to the shape of the n-string 
wave packets in real space. By inserting the complex 
momenta of the individual constituents of the string so¬ 
lutions to the Bethe wave function, it can be shown 
that for example a 2-string state with momentum 7r/2 
contains exponentials in the Bethe wave function read¬ 
ing , with an effective binding length ^ = 

2/ln(2A^). This average distance between the con¬ 
stituent particles in the string states provides a lower 
bound to how localized the wave packets constructed out 
of such string states can be. 

For a < the wave packet will start to lose its 

Gaussian shape in real space in favour of a simple ex¬ 
ponential decay around the wave packet center. This is¬ 
sue can be circumvented by choosing large enough Gaus¬ 
sian widths, but will require much higher system sizes. 
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The effective binding length becomes smaller at higher 
anisotropy, making the problem of the extra exponen¬ 
tially decaying shape to the Gaussian wave packet con¬ 
struction relevant only at |A| < 1. 

Induced by the nonlinear dispersion relation of the 
magnons, the width of a wave packet in real spa ce w ill 
furthermore increase in the course of time (see Eq. (A21), 


/ q /2 A 2 +2 

= + (28) 

obtained by expanding the dispersion relation (Eq. 0) 
to second order around the average momentum, yield¬ 
ing 6n = ^ \p=p- The broadening of n-string wave 

packets in real space is therefore described by the initial 
width a and 




( 02(0) ^ ^ 

U2n(0) 


cos^(p). 


(29) 


To first order, all strings of arbitrary length are stable 
at momenta p = ±7r/2, but possess a non-trivial de¬ 
pendence on the anisotropy for other momenta. More¬ 
over, the broadening of 1-string wave packets is not in¬ 
fluenced by the anisotropy. The anisotropy dependent 
factor can again be approximated in the large anisotropy 
limit given by Eq. (251, showing that the stability of the 


wave packets increases with increasing anisotropy and 
string length. 

In Eig. [2 the magnetization profile of a diffusing 2- 
string wave packet with zero group velocity computed 
from algebraic Bethe ansatz is shown. We used as average 
momentum p = tt since the energy dispersion relation 
has its minimum there. Furthermore, fitted parameters 
on the time-dependent wave packet width are compared 
to the theoretical values of for 2- and 3-strings. 

The diverging behaviour at low anisotropy of both 
and the effective binding length of strings constrain the 
applicability of scattering theory results for the planar 
regime |A| < 1. 


B. Stability of spin blocks 

A similar or even more drastic dispersing behaviour 
can be observed for blocks of n adjacent sites. Fig. 
shows the time evolution of a block of 20 upturned spins 
in a ferromagnetic chain of downturned spins for different 
anisotropy parameters A. The results shown in Figs. 
and were obtained using the tim e evolving block deci¬ 
mation (TEBD) algorithm .liSlIlIIll 

In the axial regime A > I, the initial state mostly 
projects onto many-magnon bound states, namely n- 
strings with n < 20, which for such anisotropy values are 
tightly bound and thus have a large overlap with the ini¬ 
tial spin block. A quantitative analysis of the overlap be¬ 
tween the initial spin block and large strings is provided 
by Refs. HHI and [501 for a comparable situation involving 



FIG. 1. Left: Decaying 2-string wave packet with zero group 
velocity (bound state of two magnons with Gaussian momen¬ 
tum distribution centered with 0 = 4 around p = tt, where 
the energy dispersion has its minimum), in a chain of A = 100 
sites, calculated from algebraic Bethe ansatz matrix elements. 
Right: prefactor 5 ^ of the d epe ndence of the wave packet 
width. Theoretical curve, Eq. (291 with p = tt and J = 1, as 
a function of A. The data points are retrieved by fitting 5 ^ 
from the decaying wave packets from the Bethe ansatz time 
evolutions. 


a prepared domain wall state containing M consecutive 
down-spins on a polarized background. The overlaps be¬ 
tween the M-spin block state and a few string configura¬ 
tions are considered, where the normalization saturation 
becomes entirely dominated by the M-string states with 
increasing anisotropy. 

Since the spin blocks mostly contain large strings 
for large anisotropy, they display slow dispersion 
(see Eq. 0 ), meaning that the initial spin block re¬ 
mains more or less intact in time over long time scales. 
However, as the isotropic point A = 1 is approached, 
the nature of the overlaps drastically changes. Eigen¬ 
states with combinations of smaller strings start carry¬ 
ing a larger fraction of the total overlap with the ini¬ 
tial state. Under time evolution, one thus sees spacetime 
propagation lines corresponding to shorter strings, which 
disperse more rapidly. By the time one has entered the 
planar regime, 0 < A < 1, the initial spin block decom¬ 
poses into all available string lengths including the most 
rapidly-dispersing I-string states, leading to a rapid dis¬ 
persion of the magnetization throughout the ‘light cone’ 
defined by the maximal group velocity of the I-strings. 

A comparison of Figs and (both at p = tt) shows 
that spin blocks decay faster than n-string wave packets 
of the previous subsection, which has a twofold explana¬ 
tion. First, an n-string wave packet is a superposition of 
n-string Bethe states which have no decay channel into 
strings of smaller lengths. In contrast, in the superpo¬ 
sition of the initial spin block state all smaller string 
lengths are allowed and actually present. Second, in 
the momentum distribution of the spin block, momenta 
belonging to high velocities are not Gaussian-like sup¬ 
pressed. Therefore, states with high velocities (p « 7r/2) 
are more dominant in the spin block state than in the 
n-string wave packet both with p = tt. However, the 
spin-block state and the n-string wave packets become 
similar at higher A, such that the analytic predictions 
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position 


FIG. 2. Decay of a block of 20 upturned spins in a chain of 
N = 1000 sites for different values of anisotropy, computed us¬ 
ing TEBD. Around the isotropic point A = 1, the behaviour 
of the spin block changes drastically, as the spin block be¬ 
comes less tightly bound and will start to decay into smaller 
strings. 


on scattering displacement becomes applicable to both 
cases. 


C. Scattering n-string wave packets 


The pre-scattering bound magnon initial states |d>(0)) 
are composed of two Gaussian wave packets, where the 
construction relies on allocation of individual momenta 
to distinct strings within a single Bethe state according to 
Eq. (14). We localize two Gaussian wave packets labeled 
by j = 1,2 with average momenta pj around two well- 
separated lattice sites 


14/(0)) = Wo ^ (pi), A(-)(P2)), (30) 

Pl,P2 

where 

^ _ g-i(PlSl+P2E2)-/^(Pl-Pl)^-/^(P2-P2)^^ 

For simplicity, we label Bethe states here only by the 
two string centers A*^"^ and A*^™^ instead of the whole 

set {A} ^ = {aS}”=i U {Ag)}r=i. The two 

centers are respectively uniquely determined by the mo¬ 
menta Pi and p 2 via Eq. ( [T^ . Due to the energy disper¬ 
sion ( [IT] ) of the bound magnons, the relative velocity of 
the wave packets is maximized at Pi = —P 2 = —7r/2. 

Fig. shows time-dependent magnetization profiles for 
scattering of 1- and 2-string wave packets respectively. 


0.5 

0.4 

0.3 

0.2 


FIG. 3. Time evolution of {Sj{t)) illustrating scattering n- 
string wave packets, computed using algebraic Bethe ansatz, 
for N = 100, A = 2, pj = —pj = —7r/2 and a = 4. Left: 
scattering of 1-string wave packets (single magnons). Right: 
scattering of 2-string wave packets (bound magnons). 



computed from the algebraic Bethe ansatz matrix ele¬ 
ments described in Sec. Eel A close examination of the 
profiles shows distinctive features akin to soliton scatter¬ 
ing, namely that the wave packets emerge out of a col¬ 
lision intact, but spatially displaced. This displacement 
will be quantified in the next section. 


D. Direct phase shift measurements 

The real-time scattering trajectories in this work are 
analyzed using the idea that the scattering between two 
localized wave packets is well-described by the scatter¬ 
ing phase shift corresponding to the average momenta 
of the two wave packets. This idea can also be veri¬ 
fied through direct observation of the phase accumulated 
by the wave function of a finite chain during a scatter¬ 
ing process. We therefore compare the evolution of an 
interacting chain (A ^ 0) with the evolution of a non¬ 
interacting chain (A = 0), each containing two localized 
1-string wave packets. The overlap between the two wave 
functions, 

rA(t) = (4/[°l(t)|4-[^l(t)), (32) 

gives the phase acquired due to the interaction between 
the magnons. The value of A is here indicated in the 
superscript. 

In Appendix]^ we show how the phase of the quantity 
Ta after a single scattering event, calculated using nu¬ 
merical exact diagonalization, matches the Bethe ansatz 
phase shift of Eq. (§, even for wave packets that are spa¬ 
tially well localized. Such overlaps between time-evolved 
wave functions are currently not directly accessible by 
Bethe ansatz. 


IV. SCATTERING DISPLACEMENT 

The initial state constructed from Bethe states is pre¬ 
pared as two wave packets of bound states of an arbitrary 
finite number of magnons with initial average positions 
























and momenta {xj, pj). The m- and n-string wave pack¬ 
ets are constructed separately at large separation, such 
that their motion before and after scattering can be con¬ 
sidered to be free. In particular, the motion of the center 
of each wave packet is (see Appendix 


. . _ 1 + '^jt before scattering, 

' ' Xj + Vjt — Xj{Pi,P 2 ) after scattering. 


(33) 

given in units of lattice distance, where the velocity Vj 
was defined in Eq. (24). Note that, in the case of single 


magnon scattering for example, a negative average mo¬ 
mentum m yields a positive velocity Vj and vice versa 
(see Eq. ([^ with J positive). 

The displacement Xi(Pi)P 2 ) can be obtained by ex¬ 
panding the scattering phas e shift around the average 
momenta, see also Eq. (|B11 ), and is therefore given as 


XjiPl,P2) = 


(34) 


Due to Eq. ^ , expression (35) for the displacement con¬ 
sists of a sum of the momentum derivative of the func¬ 
tions 9s of which the individual terms are computed as 


desjX^^Hpi) - X^^Hp 2)) _ dOsiX^^'i - X^^'i) dX<^^\pi) 

dpi dpi 

__ sin«) sin(nC) 

cos(sC) - ch(2Aj”^ - 2A^'"^) cos(nC) - cos(pi) ’ 


For the scattering phase shift of two 1-string wave 
packets withp^^ = —p 2 = —7^12, the displacement is given 
as a function of anisotropy as 


( 1 , 1 ) 

Xi 


(- 


TT TT 
2 ’ 2 / 


sin(2^) tan(^) 

cos(2^) — ch(4atanh(tan |)) 


1 + A^ 


(38) 


where we introduced the notation with subscript j to 
refer to the momentum derivative Xj = of ffi® scat¬ 


tering phase shift x- 

Eq. (33) assumes that all scatterings occur without 
particle production, which is the case for the integrable 
model we are dealing with. 


A. Displacements from Bethe ansatz 


An analytic expression for the displacement as a func¬ 
tion of anisotropy and incoming momenta can be ex¬ 
tracted from the Bethe ansatz scattering phase. The 
phase shift of two bound magnons of arbitrary length 
is obtained from the scattering kernel 0nm of the Bethe- 
Gaudin-Takahashi Eqs ([^, which consists of a sum over 
the functions 9s, s = |n —m|,... ,n-|-m, defined in Tab.|lj 
see also Eq. ([^. 

For the scattering of an n-string wave packet labeled 
by 1 with an m-string wave packet labeled by 2, the dis¬ 
placements on the trajectories of the n- and m-string 
wave packets (j = 1,2 respectively) are given as a func¬ 
tion of momenta as 


xi"’'"^(Pl,P2) = 


dOnm (A(")(pi)-AM(p2)) 


dpj 


P1=P1 

P2=P2 


(35) 

The n- and m-strings with centers A^"^ and A^™^ carry 
momenta Pi and P 2 , respectively. The wave packets are 
located such that xi <C X 2 and vi > V 2 - 

We first discuss the planar case |A| = |cos(C)| < 1. 
By inverting Eq. (12), we express the string center of an 


n-string in terms of its momentum as 
A*-"^ (p) = atanh ^tan ^ tan 


TT — p 


(36) 


The latter result could have been also obtained directly 
by taking the derivative of Eq. (§. 

Similarly, the scattering displacement of two 2-string 
wave packets becomes 


Xi 


( 2 . 2 ) 


/_ ^ _ 

V 2 ’ 2 / “ 


2sin(2C) tan(2^) 


cos(2^) — ch(4atanh(tanC)) 
sin(4(() tan(2C) 


-f 


cos(4C) — ch(4atanh(tan^)) 
16A® -4:A^ + 5 


(4A4 -K 1)(1 -h (2A2 - 1)2) 


-3. (39) 


The validity of the latter equation only extends to the 
region where A > \j\j2, as 2-strings with mom entum 
p = ±7r/2 do not exist for lower anisotropj^i^I^, which 
can be shown from the anisotropy dependent maximum 
string quantum numbers. 


For the regime A > 1, Eqs (38) and (39) hold as well. 


since the 0s (A) for both regimes are just rotated in the 
complex rapidity plane with respect to each other. Start¬ 
ing from 0s (A) for A > 1 with C, = acosh(A) therefore 
yields identical results for the scattering displacements. 

The displacements for the scattering of an 1-string 
wave packet at a 3-string wave packet is given by 


(1,3) 
Xi 


/_ _ _ 

V 2’ 2/ ~ 


(1,3) 

X2 


(1,3) 

Xl 


'_^ _ 

2 ’ 2 / ~ 


(40) 


(41) 




(4A2-1)2(4A2-3) 
2(4A2-hl)(4A4-3A2-i-l) ’ 

(4A2 - 1)3 

2(4A2 -h 1)(4A4 - 3A2 -h 1) ’ 
2A(2A-)-l)3(2A2-hA-l) 

' (2A2 -h 2A -H 1)(8A4 -h 8A3 -h 1) ’ 

(42) 

2A^{2A + 1)* 

(2A2 -f 2A -f 1)(8A4 -h 8A3 -h 1) ' 

(43) 
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B. Comparison of scattering theory and time 
evolution 


The displacement in the trajectories induced by scat¬ 
tering effects is easily deduced from the time evolution 
data of Fig. For symmetric cases with identical parti¬ 
cles, the average position of the magnetization of a single 
wave packet can be computed on one half of the system 
as a function of time, 


OXf (0 = 




Efii (i - (sm) 

^ JV/2 

i=i 




(44) 


The result for the scattering between two 2-string wave 
packets is plotted in Fig. Note that the plotted av¬ 
erage location of the wave packet in Fig. by means of 
Eq. (44) does not resemble the actual trajectories when 


the wave packets spatially overlap with each other, as 
the trajectories are not properly defined during the scat¬ 
tering event. A linear fit with the same slope is applied 
to the propagation of the average position of the wave 
packet before and after scattering. The horizontal differ¬ 
ence between the two straight lines is the displacement of 
the wave packet due to scattering effects. The procedure 
of measuring the displacements by means of Eq. (44) was 


performed for multiple values of anisotropy A. Moreover, 
these results can be compared with Eq. ( [^ , where the 
displacements for the situations of Fig. |^are specifically 
given by Eqs (38) and (l3^ as function of anisotropy. The 


results are plotted in Fig. as well, showing agreement 
between both approaches. The explicit time evolution 
relying on algebraic Bethe ansatz matrix elements pro¬ 
vides confirmation of the analytical predictions for the 
displacements. 

Besides symmetric scattering situations, colliding dis¬ 
tinct n-string and m-string wave packets can be con¬ 
structed and traced in the time-evolved magnetization 
profile. Fig. shows the scattering between an 1- and a 
3-string wave packet. Two situations are distinguished, 
the former with both wave packets at maximal veloc¬ 
ity at = —p 2 = — 7r/2, where the larger string moves 
much slower because of its effective mass, see Eq. (25). 


The latter situation consists of an incoming 1-string wave 
packet scattering on a stationary wall of a 3-string wave 
packet at |j 2 = tt. The corresponding analytic scattering 
displacements, Eqs (40)-([4^, are shown adjacent to the 
time evolution plots in Fig. The scattering displace¬ 
ments are measured from the time evolution by impos¬ 
ing a Gaussian fit on the wave packet after scattering and 
comparing the average location to the time evolution of a 
single wave packet without scattering. In the lower right 
panel (j >2 = tt), the fitting procedure of the average lo¬ 
cation of the 3-string wave packet becomes less accurate 
for decreasing A. 



Single magnons + 
Bound magnons x 



FIG. 4. Left: Measurement of the scattering displacement 
from Bethe ansatz time evolution of bound magnons at A = 
2, by computing the average location of the wave packets 
according to Eq. (44 1 and taking the horizontal difference of 
the linear fits. Right: Scattering displacement (in units of 
lattice distance) for single and bound magnons as function of 
anisotropy. Measured data from algebraic Bethe ansatz time 
evolution of magnetization (see left panel) compared to the 
derivative of scattering phase shifts, see Eqs (38)-(39). 


In the planar regime where |A| < 1, we encounter sub¬ 
stantial limitations (as described in Sec. IIIB) on both 


the construction of the scattering wave packets, as well as 
on the comparison to results of scattering theory. Due to 
the effective binding length of the individual constituents 
of the string states, the tails of the magnetization profile 
of the wave packets start overlapping with each other sig¬ 
nificantly at low A, invalidating important assumptions 
of scattering theory which include asymptotic separation 
of the wave packets before and after scattering. Compari¬ 
son of the measured displacements to scattering theory is 
therefore not meaningful for higher strings in the planar 
regime. Only going to much larger system sizes would 
resolve the aforementioned issue. 


C. Displacements from TEBD calculations for 
scattering of spin blocks compared to the Ising limit 

The scattering displacement turns out to have a par¬ 
ticularly simple form at large anisotropy. In particular, it 
was found in Ref. EH that when a propagating n-particle 
cluster is incident on a larger block, the block is displaced 
by 2n sites. This can be explained from the Bethe ansatz 
results presented in previous subsections, by taking the 
Ising limit A —>■ (X) of Eq. (35) for the displacement of 


an n-string wave packet scattering at an m-string wave 
packet. 

First, we obtain for all s 


90«(A(")(pi)-AW(p2)) 

lim -^ - ^ = — 1. 

A—yoc Opi 


(45) 


which is independent of n, m, pi, and p 2 . Using Eq. ([^ 
for the phase shift between an n-string and an m-string 
eventually yields 

lim Xi^’"‘\pi,P 2 ) = -2min(n,m) -b Snm ■ (46) 

A—>-oo 
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FIG. 5. Left: time evolution from algebraic Bethe ansatz 
of {Sj{t)} for a single magnon wave packet with momentum 
Pi = —7r/2 scattering against a 3-string wave packet with 
momentum pj = 7r/2 (top) and = n (bottom) respectively, 
with N = 100, A = 2 and a = 4. Right: Corresponding scat¬ 
tering displacements (i n un it s of lattice distance) as function 
of anisotropy, see Eqs (40 1-(|43|, where the points are mea¬ 


sured displacements obtained from the time evolution data. 


Thus, the scattering displacement for unequal wave pack¬ 
ets (n m) is equal to twice the number of particles in 
the smaller wave packet. 


The leading order term of the Ising limit in Eq. (46) 
can be given for the case where cos(p^ —p 2 )^ 0 as 


(n.m) 

Xl,LO 


— ^nm ( (2 ()nl) 


COs(Pi -P 2 ) _ . -2 

2 2 J 

X \ P2) A —In—ml 

G Onm) 2|n-m|-l ^ 

- (2 min(n, m) - 6nm) ^ , (47) 


yielding an error estimate for Eq. (46). A systematic 


expansion for all momenta becomes cumbersome due to 
the summation in Eq. 0 and is left out of consideration. 
If cos{pi —P 2 ) = fhe leading contribution is given by 
the third line of (47). If further cos(p 2 ^) = 0, the leading 
order terms will be formed by higher powers like 
for n = m or for m ^ n. 

Fig.|^ shows TEBD results of the scattering of a two- 
spin block excitation with a block consisting of 10 adja¬ 
cent spins at A = 5.0. The initial 2-spin block was cre¬ 
ated by upturning two neighboring spins at lattice sites 
2 and 3. The first site (with open boundary conditions) 
is energetically inaccessible for large A; hence the 2-spin 
block travels to the right and is incident on the 10-spin 
block. The displacement is clearly by 4 sites, as predicted 
by Eq. (46). A similar displacement by 2 sites in the case 


of a single incident particle was highlighted in Ref. Uni 



10 20 30 40 50 

position 


FIG. 6. TEBD time evolution of scattering of a two-spin block 
state with a block of 10 upturned spins at A = 5. The block of 
10 upturned spins is shifted by four lattice sites upon impact 
of the 2-string like bound magnon. This is in correspondence 
with Eq. (461, yielding a shift of four caused by a 2-string on 
a larger string. 


Although the 2-spin block and the 10-spin block are 
not explicity prepared as wave packets in this case, at 
large A, string states are tightly bound, and therefore 
these blocks may be understood intuitively to be close to 
2-string wave packets and 10-string wave packets. The 
Ising limit Eq. (46) thus provides a satisfactory explana¬ 
tion to the shift observed in these numerical experiments. 

Even more, it was observecP^ that the behaviour of 
the scattering displacement of a large spin-block state 
close to the isotropic point at A = 1.1 still resembles the 
scattering behaviour of the large A limit by shifting the 
spin-block by two lattice sites upon scattering. In order 
to explain this from Bethe ansatz, we take the limit of 
Eq. (35) for the scattering of a large string (m S> 1) with 
center A^"*^(p 2 ) and an 1-string with center A(^)(pi) at 
all values of A > 1, 


dOs (AW(pi)-A("^)(p2)) 
dp2 


m»l 


(1 + 0(6-'*'^)) (l + 0(e-'"«)) , 


(48) 


where s can only take on the values s = m—1 or s = 
m -I- 1, due to Eq. The latter equation finally gives 
for the displacement of the large m-string at A > 1, 

= 2 + 0(e-(™-i)«), (49) 

implying that already at A = 1.1 = 0.4436), the scat¬ 

tering displacement of a large spin block should still be 
close to two sites, as is in agreement with the aforemen¬ 
tioned observation. 


V. CONCLUSIONS 

In this work, we have studied the quantum analogue of 
soliton-like scattering phenomena in the anisotropic spin- 
1/2 Heisenberg chain, by utilising the algebraic Bethe 
ansatz. We considered quantum scattering of localized 
excitations, created from linear combinations of Bethe 
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states with Gaussian-distributed momenta, constructing 
wave packets of n bound magnons. This construction 
allows to study scattering phenomena of wave packets 
containing an arbitrary number of bound magnons. 

Exact methods based on the algebraic Bethe ansatz 
provide a framework to evaluate the time-dependent ex¬ 
pectation value of the local magnetization profile, which 
allows for a spatial tracking of the localized excitations. 
This explicit unitary time evolution of the initial state 
relies on the availability of determinant expressions for 
matrix elements of local spin operators. 

The algebraic Bethe ansatz time evolution of colliding 
wave packets of bound magnons displays a spatial dis¬ 
placement in the trajectories of the wave packets under 
scattering, consistent with scattering theory results. For 
different values of anisotropy, fits on the displacements of 
the time evolved trajectories are in agreement with an¬ 
alytical results on the displacement from the derivative 
of the Bethe ansatz scattering phase shifts, for several 
combinations of string lengths. 

The scattering phase shift can also be measured di¬ 
rectly as well for the scattering between two localized 
single-magnon wave packets, again matching phase shift 
expressions provided by Bethe ansatz. Using TEBD, 
scattering displacements from spin-block states have 
been studied, showing similar scattering features and val¬ 
idating the analytic predictions of the Ising limit for the 
scattering displacement. 

The experimental realizability of real time tra cking 
of localized excitations in the Heisenberg spin chairP^^l^i] 
might provide an opportunity to study dynamical scat¬ 
tering phenomena of (bound) magnons. A possible man¬ 
ifestation of such phenomena might be provided by the 
soliton-like scattering effects analysed in this work. 

The results on the scattering displacements can be ex¬ 
tended to other Bethe ansatz solvable models. Finally, 
the time evolution method relying on matrix element ex¬ 
pressions from algebraic Bethe ansatz can be used to con¬ 
struct other initial states in spin chains as well and to 
study their respective relaxation phenomena. 
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Appendix A: General results from scattering theory 


In this Appendix, we review some general results 
from quantum scattering theory, emphasizing the direct 
connection between physically observable quan tities and 
(derivatives of) the scattering phase shift.^IEll 

We consider an initial state with two well-separated 
quasiparticles (e.g. magnons and magnon bound states) 
with almost well-defined positions (aii, X 2 ) and momenta 
(Pi, P 2 )- In III® asymptotic region, i.e. when the distance 
between the two wave packets is much larger than the ra¬ 
dius of the interaction, the time evolution is free. Hence, 
the centers of the wave packets Xj{t) = {xj{t)) translate 
rigidly. The scattering between the two particles however 
introduces a displacement proportional to the derivative 
of the scattering phase x(PijP 2 )- More precisely, in the 
asymptotic regions, the motion of the center of each wave 
packet is 


. X _ J + Vjt before scattering, 

' \ Icj + Vjt — Xj{PiTP 2 ) after scattering, 

(Al) 

where Xj = ^ is the displacement, while the group 
velocity Vj is given by the derivative of the dispersion 
relation, i.e. Vj = 

Similarly, the scattering has an effect on the width of 
each wave packet ^x‘j{t) = {[xj{t) — Xj(t)]'^). For Gaus¬ 
sian wave packets we have to first order 


Axf{t) = 


QU I 

4 + a? 


X12 


(xii-t<5i) 

-- 


before scatt., 
after scatt., 


while Sj = 


dEj 


(A2) 
_ and 

t =0 ' 'P3-P3 

An analogous formula holds 


where aj = (f) 

y = I 

dpidpj \pi=p.,pj=p^ 

for Ax 2 {t) as well. 

Scattering also builds up correlations between the 
(initially uncorrelated) Gaussian wave packets, as can 
be seen from the time evolution of the correlator 

AxiAx2{t) = ([il(t) - Xi{t)] [X2{t) - X2{t)]), 
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AxiAx2it) 


f 

0 

[ X 12 

Xll—tSl 1 X22 — tS2 

-^5- “T - ^ - 

[ af ^2 J 


before scatt., 
after scatt. 


(A3) 


All the aforementioned quantities carry information 
about the derivatives of the scattering phase and can be 
in principle measured in a scattering experiment. These 
results hold for any one dimensional theory with stable 
particles. For the XXZ spin chain the stability of magnon 
bound states above the binding energy threshold is pre¬ 
served by the integrability of the theory. 
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Appendix B: Scattering of two particles in one 
dimension 

In this Appendix, we review some general results for 
the scatteri ng of two particles (magnons, magnon bound 
states, etc.),*^i^ and derive Eqs ( |Al| )-( [X^ . For simplic¬ 
ity, we consider the scattering of two distinguishable par¬ 
ticles in a continuum integrable model. The same results 
can be obtained for identical particles. 

So, let us consider two particles with asymptotic mo¬ 
menta pi andp 2 and different dispersion relations Ei{pi), 
i = 1,2 in an infinite volume (zero density). The state¬ 
ment of the coordinate Bethe ansatz, see Eq. § , is 
that in the asymptotic region where the two particles 
are very far apart the eigenfunctions of the system are 
plane waves. 


, xi<X2, 

j »X2, 

(Bl) 

where S{pi,p 2 ) = —e^^^'Pi,P 2 ) jg scattering matrix and 
x(pi,P 2 ) the scattering phase shift. At zero density, the 
energy is simply 

E{pi,P2) = Ei{pi) + E2 {p2) ■ (B2) 


Let us briefly comment on the structure of the wave func¬ 
tion (Bl). First of all, as we discussed in the main text 


for the XXZ model, bound states are characterized by 
complex conjugate rapidities, which leads to exponen¬ 
tially decaying terms in the Bethe wave function (§ with 
respect to the relative coordinate. This is a feature of 
bound states that are characterized by a center of mass 
coordinate. In what follows, we do not denote these ex¬ 
ponentially decaying terms and label the bound states 
only with the position of the center of mass. For elemen¬ 
tary particles, Eq. (Bl) is a consequence of the conser¬ 


vation of energy and momentum in one dimension, and 
as such it is valid for any model with a sufficiently short- 
range potential. Instead, if one of these particles is not 
elementary but it is a bound state, then for a general 
theory the previous simple form of the wave function is 
not true anymore. The bound state can decay and scat¬ 
tering be diffractive. However, there exist models (such 
as the XXZ spin chain) for which the scattering is always 
non diffractive. Hence, the coordinate Bethe ansatz de¬ 
scribes a complete set of asymptotic eigenfunctions, as 
thoughtfully discussed in Sutherland’s book.E^l For such 
theories, bound states cannot decay, but are protected 
by integrability. 

Let us consider the scattering problem. At time t = 0 
the two particles are far apart and have (almost) well- 
defined positions Xj and momenta j = 1,2. Without 
loss of generality, we may assume that Xi <C X 2 and vi > 
U 2 , where vj is the group velocity. 


dE, 

dp. 


(B3) 


Pj=Pj 


When they are far apart, the two wave packets move 
with velocities vj. If vi < V 2 , the evolution is always 
free, while for vi > V 2 at some time the two particles 
become close and the interaction plays a role. 

The time evolution of a two-body wave function is 
given by 


tp{Xi,X2,t) 


where 

C{pi,P2) 


f ^ ^ 

J 27r 27r 


C{P1,P2) 


^-it{E{px)+E{p2)) 


‘‘Ppi ,P2 (^1 5 ^ 2 ) 


(B4) 


J dxidx2Pp,^p^{xi,X2)'<p{xi,X2), (B5) 


Pp,,P 2 is a complete set of eigenfunctions, and 'ip{xi,X 2 ) 
is the initial wave function. Now, we consider an initial 
state which is factorized. 


tp{xi,X2) = 'tpl{xi)lp2{x2) (B6) 

where 

= 1^ Up,) , j = 1 , 2 . (B7) 


Since we are interested in computing also the spreading 
in time of the wave packet, we assume that the functions 
il)j{xj) are Gaussian wave packets with minimal indeter¬ 
mination, 

EP^^^e^ , (B8) 


where Aa;| = f dxj (x, — Xj)^ li/^j(xj)l^ = <a|/4. This 
assumption may be dropped if we are interested only in 
the displacement, Eq. ( |A1[ ). 

Up to now, everything is exact. However, we can make 
the following approximations: 


• Since the initial state ip{xi,X 2 ) is sharply peaked 
around xi and X 2 w ith xi <tiX 2 , the only relevant 
contributions to (B5) come from the regions around 
these two poi nts. We can then use the asymp¬ 
totic formula (Bl) for (ppj^(a:i, 2 : 2 ) with xi X 2 - 
Inserting Eq. (Bl) into (B5), and using definition 
(B7) yields 


C{pi,P2) « Uvi)Up2) 


(B9) 


Since the functions 'ijjj{pj), j = 1,2, in ( |B9| ) are 
peaked around pj « Pj , we can expand the disper¬ 
sion relations Ej( pj) a round these momenta in the 
integrand of Eq. ( |B4[ ). Moreover, for xi ^ X 2 we 
may substit ute the asymptotic expression (Bl) for 
Ppi,P 2 iiito (B4) and similarly expand the scatter¬ 


ing phase shift xiPijP^)- Up to second order these 
expansions read 


EjiPj) = Ej{p,) + Vj{p, - pA + -^{pj - p,)= 


(BIO) 
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FIG. 7. Left: Time evolution of (5|(f)), performed using 
numerical exact diagonalization for L — 26 sites. Two single- 
magnon wave packets are prepared with initial width a = 
5/\/2 « 3.5355 and opposite momenta ±fc; here k = 0.57r 
and A = 2.5. Right: The Bethe ansatz phase shift (lines) 
are compared with the phase of T (symbols, shifted by tt), 
obtained from the wave function of the time-evolved L = 26 
chain. The comparison is shown for three k values. 


with Sj 


I 

dpj \pj=Pj’ 


and 


X{Pl,P2) = X{Pl,P2) + Xz{Pl,P2)iPi - Pz) 

+ \{P^-Pi)Xij{Pl,P2)iP3 -Pj) (Bll) 


with Xi = and Xij = Opidpj ’ 
repeated indices. These expansions up to the sec¬ 
ond order are physically meaningful only if we can 
ignore the distortion of the wave packet, and so 
they are no more valid for times long enough. 


Taking advantage of these approximations, we are now 
in the position to derive Eqs (Al|-(A3|. Before scatter¬ 
ing, the two wave packets propagate freely. Since they 
start around xi and X 2 with xi X 2 , they are cen¬ 
tered around Xj{t) with Xi{t) <C X 2 (t) for small times. 
Thus, we can use the asymp toti c formula of i^pj_p 2 (a:i, 0 : 2 ) 
valid for xi <C X 2 in Eq. (Bll, and hence perform the 
Gaussian integrations thus obtaining the pre-scattering 
results, Eqs (Al)-(A3). Similarly, long after the scatter¬ 
ing, we have xi{t) ^ X 2 (t). Th us, t aking advantage of 
the proper asymptotic formula (Bl) for the eigenfunc¬ 
tions and the expansion (Bill of the scattering phase 
shift we also obtain the post-scattering formulas (Al)- 
dASl). 


Appendix C: Obtaining the phase shift from the 
phase of the wave function of the full chain 

A prominent theme of this work is that real-time 
scattering between two localized wave packets is well- 
described by the phase shift corresponding to the average 
momenta of the two wave packets. In this Appendix we 
connect real-time scattering data for two single-magnon 
(1-string) wave packets in a hnite chain directly to the 


definition of the scattering phase shift, namely, that the 
scattering phase shift is the phase picked up by the sys¬ 
tem wave function during the scattering process. To 
demonstrate the robustness of this idea, we show results 
from a stringent situation of rather small wave packets 
(width of a few sites) in a rather small chain (« 25 sites), 
far from the usual idealized limit of infinitely extended 
excitations. 

Here the Hamiltonian of the XXZ model is used in the 
following form: 


(Cl) 

The extra terms compared to Eq. 0 are convenient for 
considering the phase of the time-evolving wave func¬ 
tion in the two-magnon sector. When evolving with the 
Hamiltonian ([^, there is a constant accumulation of A- 
dependent phase in the time evolution even when the 
magnons are spatially separated; this is avoided with the 
above form. With this form, the interaction affects the 
phase of the chain wave function only when the magnons 
meet each other. 

In the initial state, each magnon wave packet is pre¬ 
pared as a Gaussian, localized approximately around L/4 
and 3L/4 respectively, with opposite momenta ±/c. The 
initial state is thus 


L-l 


i/A = jE 


7 = 1 


IV'(O)) = A'oE'5'a; exp 


{x - XqY 


^ — ikx 


Y^Sy exp 


{y-yof 




^+iky 1 ^^ ^ ^ ^ 


(C2) 


with xq and 7/0 near the center of the left half and right 
half of the chain. Here x and y are used as discrete 
site indices, and Afo is a normalization constant. The 
definition of the width a is chosen to be consistent with 
a in the main text. To avoid excessive dispersion, we 
use k near tt/2. This preparation ensures that the two 
wave packets will move toward each other, collide, and 
continue on after the scattering, as shown in the left panel 
of Fig. The effect of the interaction JA should be 
felt only when the particles are crossing each other. This 
interaction gives a phase shift to the wave function during 
the scattering. 

The wave function at time t is 


|7/;Al(t)) = e-*^^‘|7/>[^l(0)) = e-*^^‘|7/i(0)) (G3) 

(same initial state for every A). We consider the overlap 
rA(T) = (7/;[°l(r)|7/;[^l(T)) (G4) 


at some appropriately chosen ‘final’ time t = T. Evolu¬ 
tion with the A = 0 (non-interacting) Hamiltonian gives 
a ‘reference’ state. We are interested in the phase accu¬ 
mulated in the evolution with the A 7 ^ 0 Hamiltonian in 
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comparison with the reference state, namely the phase of 
Fa- The final time T is chosen such that the particles 
have completed their scattering, but have not reached 
the boundaries of the chain. So, there are no edge ef¬ 
fects. For example, in the process shown in Fig. [^left), 
it would be reasonable to compare overlaps at T^13J“^. 

The phase of Fa should approximate the scattering 
phase shift. 


lim arg(FA) = % - tt . (C5) 

L—>-oo,a—>-oo 

The shift tt is in accordance with the convention used in 
this paper, e.g., in Eq. For k = k 2 = —ki, the Bethe 


ansatz phase shift is obtained from Eq. (§ to be 


y = TT — 2 atan 


/ Asin(fc) \ 
\1 — Acos(fc) / 


(C6) 


Fig. [fright) compares the phases obtained from the 
time-evolving wave function of a finite chain contain¬ 
ing two relatively narrow wave packets, with the Bethe 
ansatz phase shift expressions which are strictly valid for 
delocalized excitations. Even with wave packets as nar¬ 
row as a « 3.5, the Bethe ansatz expressions match ex¬ 
tremely well the phase acquired in real-time evolution. 
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